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I .  INTRODUCTION 


Recently,  there  has  been  a  renewed  interest  in  bombs  impacting  concrete 
targets  and  their  trajectory  through  the  target.  This  is  especially  important 
when  the  bomb  is  equipped  with  a  delay  fuze.  Experimental  firings  conducted 
by  Roecker1  of  1/lSth  scale  MARK84  (MK84)  bombs  impacting  thick  concrete  tar¬ 
gets  at  various  obliquity  angles  and  striking  velocities  have  shown  the  bomb 
to  follow  one  of  two  trajectories: 

(a)  penetration  -  piercing  of  the  target  by  the  bomb.  In  some  of  the 
tests,  partial  penetration  of  the  target  was  observed  followed  by  ejection  of 
the  bomb  from  the  crater  and  tumbling  once  clear  of  the  target.  This  behavior 
is  termed  rebound. 

(b)  ricochet  -  deflection  of  the  bomb  from  the  target  surface  in  which  the 
nose  of  the^ bomb  is  never  completely  embedded  in  the  target.  A  subset  of  ric¬ 
ochet  is  broaching  wherein  the  bomb  tunnels  through  the  target  along  a  trajectory 
parallel  to  the  target  surface,  and  after  sometime  it  turns  towards  the  target 
surface  and  exits. 

The  experiments  conducted  by  Roecker1  established  critical  obliquity  angles 
for  delineating  penetration  and  ricochet  as  a  function  of  striking  velocity  for 
solid  and  thin-walled  MK84  bombs.  The  bombs  were  inscribed  by  a  polypropolux 
sabot  and  fired  from  a  57mm  smooth-bore  launcher  at  velocities  from  150  m/s  to 
270  m/s.  The  concrete  targets  which  were  15.24cm  thick  and  approximately  1 
meter  square  were  placed  at  obliquities  from  20°  -  40°.  The  experimental  evi¬ 
dence  consists  of  post-impact  measurements,  flash  radiographs,  polaroid  photo¬ 
graphs,  and  some  movies.  Figure  1  shows  a  schematic  of  the  thin-walled  bomb. 

The  solid  bomb,  which  was  fired  to  investigate  the  influence  of  wall  thickness 
on  the  critical  obliquity  angle,  has  the  same  external  dimensions. 

Figure  2  is  a  photograph  of  typical  pre-  and  post-impact  conditions  for  a 
thin-wall  bomb  impacting  the  concrete  target  at  159  m/s  at  an  obliquity  of  31°. 
This  happens  to  be  a  ricochet  mode  and  the  ogival  nose  essentially  remains  un¬ 
deformed.  What  appears  to  be  an  indentation  in  photograph  near  the  nose  is  a 
scar  not  a  localized  plastic  deformation.  The  deformation  near  the  rear  of  the 
bomb  is  the  result  of  the  tail  of  the  bomb  "slapping"  the  target  in  the  late 
stages  of  the  ricochet  event. 

To  gain  further  insight  into  characterizing  the  mechanisms  of  the  impact 
event  which  distinguish  the  different  modes  of  penetration,  numerical  simulation 
of  scaled  MK84  bomb  impacts  have  been  performed  using  the  EPIC-32  code.  The 
results  of  these  simulations  are  discussed  and  compared  with  experimental  ob¬ 
servations. 


E.  T.  Roeoker ,  et.  al.,  BRL  Report  in  progress. 

G.  R.  Johnson ,  D.  J.  Vavriok ,  and  D.  D.  Colby,  "Further  Development  of 
EPIC-3  for  Anisotropy,  Sliding  Surfaces,  Plotting  and  Material  Models," 
Ballistic  Research  Laboratory ,  ARBRL-CR-00429 ,  May  1980  (AD  B048305L) . 
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NOTE:  ALL  DIMENSIONS  IN  cm. 

CONSTANT  WALL  THICKNESS 


Figure  1.  Scaled  MARK84  Thin 


II.  NUMERICAL  MODEL 


The  thin-walled  bomb  has  been  modelled  with  four  rings  of  elements  with  the 
outer  ring  of  elements  representing  the  2mm  steel  casing.  The  inner  rings  of 
elements  represent  the  inert  simulant  to  represent  the  explosive  tritonal.  The 
bomb  contains  855  nodes  and  3456  elements  and  the  entire  model  is  comprised  of 
8070  nodes  and  39,744  constant  strain  tetrahedral  elements. 

A  deep  penetration  "tunneling"  model  has  been  employed  for  treatment  of  the 
sliding  surface.  The  bomb  has  been  declared  the  master  and  the  target  has  been 
declared  the  slave.  This  is  the  reverse  of  that  employed  by  Johnson.* 

The  constitutive  model  in  EPIC-3  employs  an  incremental  elastic-plastic 
relationship  using  the  von-Mises  yield  criterion  to  describe  the  deviatoric 
behavior  of  materials.  Provision  is  made  to  account  for  strain  rate,  temper¬ 
ature,  and  strain  hardening  effects.  Thermal  softening  and  strain  rate  effects 
were  not  taken  into  account  due  to  lack  of  data  to  support  such  models.  The 
hydrostatic  component  of  stress  is  obtained  from  the  Mie-Gruneisen  equation  of 
state  for  the  steel  casing  and  explosive  simulant. 

A  plain  (unreinforced)  concrete  equation  of  state  which  is  valid  up  to  30 
GPa  was  used  to  model  the  concrete.  The  concrete  model  is  discussed  in  Refer¬ 
ences  2  and  4  and  is  presented  briefly  here  for  completeness.  The  concrete 
model  is  constructed  from  Hugoniot  data  by  Gregson5  and  static  yield  strength 
data  from  Chinn  and  Zimmerman.6  The  aggregate  used  in  Gregson* s  work  had  an 
average  diameter  of  3.2mm  and  an  average  initial  density  of  2.185  g/cm3,  The 
concrete  consisted  of  18%  voids,  25%  granite,  20-25%  quartz  grains,  and  25-30% 
cement  paste  by  volume. 

The  pressure-volume  relationship  for  the  concrete  model  is  shown  in  Figure  3. 
The  estimated  unloading  path  is  one  which  permits  the  concrete  to  return  to  its 
solid  density  (2.6  g/cm3)  when  loaded  past  the  point  where  all  voids  have  col¬ 
lapsed,  u  =  0.3. 


G .  R.  Johnson,  Honeywell  Inc.,  Defense  Systems  Division,  Private  communication. 

4  J.  J._  Osborn  and  D.  A.  Matusha,  "Dynamic  Response  of  a  Kinetic  Energy  Pane- 
trator ,  Vol  II.  -  Hydrocode  Analysis,"  A  FAIL  TR  78-24 ,  March  1978. 

5  V.  G.  Gregson,  Jr.,  "A  Shock  Wave  Study  of  Fondu-Frye  WA-1  and  a  Concrete," 

DNA  report  2797F,  February  1972. 

ft 

J.  Chinn  and  E.  M.  Zimmerman,  " Behavior  of  Plain  Concrete  Under  Various  High 
Triaxial  Compression  Loading  Conditions ,"  AFWL  TR  64-163,  August  1966. 
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Figure  3.  Pressure-Volume  Relationship  for  Concrete  Model 


PRESSURE  (MPa) 


The  unloading  path  is  a  straight  line  which  connects  with  the  load  curve  at  the 
point  where  all  the  voids  have  collapsed  and  has  a  slope  equal  to  that  of  the 
load  curve  at  the  point  where  total  void  collapse  occurs.  For  values  of  excess 
compression  below  that  where  void  collapse  is  initiated  (y  =  .002),  loading, 
unloading,  and  reloading  occur  along  the  same  elastic  path.  For  values  of  excess 
compression  between  these  extremes,  the  bulk  modulus  is  assumed  to  vary  linearly 
between  KQ  and 

The  concrete  model  has  been  extensively  and  successfully  employed  by  the  Air 
Force  Armament  Laboratory  for  penetrations  in  the  60  to  300  m/s  velocity  range. 4 
A  detailed  description  of  the  model  can  be  found  in  References  2  and  4.  Other 
pertinent  literature  is  found  in  References  7  and  8.  The  use  of  the  model  should 
be  limited  to  calculations  with  an  impedance  mismatch,  such  as  steel  and  concrete 
in  this  case,  because  the  model  assumes  no  transmission  of  shock  waves  from  the 
concrete  into  the  penetrator. 

In  the  calculations  presented  here  the  bomb  elements  were  not  permitted  to 
fail.  However,  the  concrete  elements  were  completely  failed,  at  which  point  the 
element  no  longer  supports  a  pressure  or  tensile  stress,  when  its  equivalent 
plastic  strain  exceeded  30  percent.  It  should  be  noted  that  preliminary 
calculations  which  raised  the  equivalent  plastic  strain  from  10-30  percent 
showed  no  distinct  differences  in  target  crater  profile. 

III.  DISCUSSION  OF  RESULTS 

Computations  were  performed  for  the  thin-walled  structure  (Figure  1)  at  obliq¬ 
uity  angles  of  25.7°,  30°,  and  33.8°  with  a  striking  velocity  of  228  m/s.  These 
obliquity  angles  and  striking  velocity  were  chosen  due  to  the  availability  of 
experimental  data.  0The  experimental  data  indicated  that  with  a  striking  velocity 
of  228  m/s  the  25.7°  impact  has  a  50%  probability  of  ricochet  and  a  50%  proba¬ 
bility  of  rebound.  The  higher  obliquity  angles  resulted  in  ricochet  or  broaching. 

It  should  be  noted  that  in  light  of  the  excessive  amount  of  computer  time 
required  by  these  calculations  the  solid  bomb  has  not  been  extensively  pursued 
and  the  results  herein  pertain  to  the  thin-walled  structure  only.  The  explicit 
integration  scheme  employed  in  EPIC-3  required  approximately  85  cp  seconds  per 
microsecond  of  simulation  on  a  CDC  Cyber  76  for  the  results  presented  herein. 

Figure  4  shows  the  deformation  profiles  for  the  three  different  obliquity 
angles  at.  the  time  where  maximum  penetration  is  achieved.  Appendix  A  contains 
a  complete  temporal  sequence  of  deformation  profiles  for  the  different  obliq¬ 
uities.  In  general,  the  deformation  profiles  in  Figure  4  closely  resemble  those 
observed  in  the  experiments.  The  discontinuity  in  the  nose  of  the  bomb  is  an 
artifact  of  the  sliding  surface  technique.  The  target  crater  profiles,  in  the 
symmetry  plane,  are  delineated  by  elements  which  have  not  exceeded  an  equivalent 
Dlastic  strain  of  30  percent. 


D.  McHenry  and  J.  Kami,  "Strength  of  Concrete  Under  Combined  Tensile  and 
Compressive  Stress, "  American  Concrete  Institute,  page  54,  10  April  1958. 

K.  E.  Crawford,  C.  J.  Higgins,  and  E.  H.  Bultmarm,  "The  Air  Force  Manual  for 
Design  and  Analysis  of  Hardened  Structures,"  AFUL  TR-74-102,  October  1074, 
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The  net  velocity  and  angular  velocity  histories  for  the  different  obliquities 
are  compared  in  Figures  S  and  6,  respectively.  As  can  be  seen  in  Figure  5,  at 
approximately  0.3ms  the  net  velocity  of  the  25.7°  impact  starts  to  increase. 
Furthermore,  the  vertical  comment  of  velocity  changes  sign.  This,  in  our  opin¬ 
ion,  is  indicative  of  a  rebound  mode.  One  can  also  note  the  slower  deceleration 
of  the  higher  obliquity  impacts. 

The  comparison  of  angular  velocity,  i.e.  the  rate  at  which  the  tail  of  the 
bomb  is  falling  towards  the  target,  histories  for  the  different  obliquities  pre¬ 
sented  in  Figure  6  show  that  the  higher  obliquity  impacts  result  inQa  large 
angular  acceleration  being  imparted  to  the  bomb.  In  fact,  the  33.8  impact  has 
been  carried  out  beyond  one  Billisecond  at  which  time  there  appears  to  be  a 
steady  state  angular  velocity  of  S00  rad/s  and  the  net  velocity  has  been  reduced 
to  approximately  10%  of  the  striking  velocity.  This  angular  velocity  history 
for  the  33.8°  impact  coupled  with  the  deformation  profile  at  1.1ns  (see  Figure 
4  which  indicates  the  initiation  of  tunneling)  indicates  a  broaching  mode,  in 
our  opinion.  The  intermediate  obliquity  is  felt  to  be  a  ricochet  mode  in  that 
no  tunneling  behavior  has  been  observed  out  to  600  ms. 

The  response  of  the  thin-walled  bomb  at  the  various  obliquities  is  further 
compared  in  Figures  7-9  which  show  the  time  history  of  bending  moment ,  shear 
stress,  and  axial  load,  respectively,  for  the  indicated  obliquities.  These 
temporal  histories  are  for  the  layer  of  elements  in  the  cylindrical  portion  of 
the  bomb  just  above  the  ogival  nose.  No  distinct  difference  is  noted  for  dis¬ 
tinguishing  penetration  and  ricochet  based  on  the  temporal  histories  for  bend¬ 
ing  moment  and  shear  stress.  However,  the  reduction  in  axial  load  beyond 
0.3ms  for  the  25.7°  obliquity  (Figure  9),  in  conjunction  with  the  correspond¬ 
ing  net  velocity  history  (Figure  6),  is  indicative  of  "rebounding." 

IV.  SUMMARY 

Shortly  after  impact,  an  angular  acceleratic.i  is  imparted  t*'  the  .>omb  and 
the  magnitude  and  duration  of  this  angular  acceleration  is  proport i  ,nal  to  the 
incidence  angle.  Eventually,  the  angular  acceleration  imparted  to  the  bomb 
reaches  a  steady  state.  A  large  steady  state  angular  velocity  appears  to  be 
characteristic  of  a  broaching  mode  of  penetration  accompanied  by  slower  deceler¬ 
ation  of  the  bomb. 

The  computer  simulations  presented  here  predict  the  critical  angle  for  rico¬ 
cheting  within  3-5  .  They  are  quite  expensive  (80-90  CPU  seconds/ 1  us  of  real 
time).  This  is  due  principally  to  the  fact  that  explicit  integration  is  used 
throughout.  While  accuracy  considerations  require  a  time  step  characteristic 
of  explicit  integration  schemes  during  the  early  stages  of  impact  (0-100  us), 
low  impact  velocity  computations  such  as  these  could  be  conducted  far  more  ec¬ 
onomically  by  switching  to  an  implicit  scheme  once  details  of  wave  motion  are 
no  longer  important.  The  principal  advantage  of  codes  such  as  EPIC-3  for  these 
applications  is  their  ability  to  supplement  the  meager  information  obtained  from 
impact  experiments.  The  judicious  combination  of  computations  and  experiments 
can  lead  to  effective  designs  with  cost  and  efficiency  factors  below  those  from 
either  technique  alone. 
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Figure  Al.  Deformation  Profile  for  25.7°  Obliquity  at  T  =  10  ps 
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Figure  A2.  Deformation  Profile  for  25.7°  Obliquity  at  T  =  20  ps 
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Figure  A3.  Deformation  Profile  for  25.7°  Obliquity  at  T  =  48.2  ys 
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Figure  A6.  Deformation  Profile  for  30°  Obliquity  at  T  =  40  us 
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Figure  A7.  Deformation  Profile  for  30°  Obliquity  at  T  =  60  us 
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Figure  A10.  Deformation  Profile  for  33.8°  Obliquity  at  T  =  40  us 
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Figure  A12.  Deformation  Profile  for  33.8°  Obliquity  at  r  =  80  ps 


-J0.00  -gO.OO  -15.00  -.10  00  -S.00  0.00  S..00  10.00  Ig.QQ  2p.ac 


YMIN  .0.00000000 


33.8  DEGREE  08LI0UITY 


TIME 

YMAX 


<.00100001 

.0.00000000 


-W.ao  -ho.  oo  -is. oo  -Woo  Too  o!ac  Too  tb.ao  Too  fb.ao  i.oo 

X  (CM) 

Figure  A13.  Deformation  Profile  for  33.8°  Obliquity  at  T  =  100  us 
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APPENDIX  B 

Input  Data  for  EPIC-3  Preprocessor  for  33.8°  Obliquity 
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2  Orlando  Technology,  Inc. 

P0  Box  855 

ATTN:  Mr.  J.  Osborn 
Mr.  D.  Matuska 
Shallmar ,  FL  32579 

1  Lawrence  Livermore  Laboratory 
ATTN:  Dr.  M.  L.  Wilkins 
P.0.  Box  808 
Livermore,  CA  94550 

1  Rockwell  International 
Missile  Systems  Division 
ATTN:  A.  R.  Glaser 
4300  E.  Fifth  Avenue 
Columbus,  OH  43216 

3  Schumberger  Well  Services 
Perforating  Center 
ATTN:  J.  E.  Brooks 

J .  Brookman 
Dr.  C.  Aseltlne 
PO  Box  A 

Rosharon,  TX  77543 

1  Science  Applications,  Inc. 

101  Continental  Building 
Suite  310 

El  Segundo,  CA  90245 

1  Lawrence  Livermore  Laboratory 
ATTN:  Dr.  J.  0.  Hallquiat 
P.0.  Box  808 
Livermore,  CA  94550 

1  Systems,  Science  and  Software 
P0  Box  1620 

ATTN:  Dr.  R.  Sedgwick 
La  Jolla,  CA  92037 
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2  TRW 

One  Space  Park,  Rl/2120 
ATTN:  D.  Auehertaan 
M.  Bronstein 

Redondo  Beach,  CA  90278 

1  United  Technologies 
Research  Center 
438  Weir  Street 
ATTN:  P.  R.  Fitzpatrick 
Glastonbury,  CT  06033 

1  US  Steel  Corporation 
Research  Center 
125  Jamison  Center 
Monroeville,  PA  15146 

1  VPI  &  SU 

106C  Norris  Hall 
ATTN:  Dr.  M.  P.  Kamat 
Blacksburg,  VA  24061 

2  Vought  Corporation 
PO  Box  225907 
ATTN:  Dr.  G.  Hough 

Dr.  Paul  M.  Kenner 
Dallas,  TX  75265 

1  Wes tinghouse ,  Inc. 

P0  Box  79 

ATTN:  J.  Y.  Fan 

W.  Mifflin,  PA  15122 

1  Drexel  Institute  of  Technology 
Department  of  Mechanical  Engr. 
ATTN:  Dr.  P.  C.  Chou 
32d  and  Chestnut  Streets 
Philadelphia,  PA  19104 

5  Southwest  Research  Institute 
Dept,  of  Mechanical  Sciences 
ATTN:  Dr.  U.  Lindholo 
Dr.  W.  Baker 
Dr.  R.  White 
Dr.  M.  F.  Kanninen 
Dr.  C.  Anderson 
8500  Culebra  Road 
San  Antonio,  TX  78228 


4  SRI  International 

333  Ravenswood  Avenue 
ATTN:  Dr.  L.  Seaman 
Dr.  L.  Curran 
Dr.  D.  Shockey 
Dr.  A.  L.  Florence 
Menlo  Park,  CA  94025 

2  University  of  Arizona 

Civil  Engineering  Department 
ATTN:  Dr.  D.  A.  DaDeppo 
Dr.  R.  Richard 
Tucson,  AZ  85721 

1  University  of  Arizona 
School  of  Engineering 
ATTN:  Dean  R.  Gallagher 
Tucson,  AZ  85721 

1  University  of  California 
ATTN:  Dr.  N.  Ziv 
504  Hilgard  Avenue 
Los  Angeles,  CA  90024 

1  University  of  California 
Department  of  Physica 
ATTN:  Dr.  Harold  Lewis 
Santa  Barbara,  CA  93106 

2  University  of  California 
College  of  Engineering 
ATTN:  Prof.  W.  Goldsmith 

Dr.  A.  G.  Evans 
Berkeley,  CA  94720 

2  University  of  Delaware 
Department  of  Mechanical 
Engineering 
ATTN:  Prof.  J.  Vinson 
Prof,  B.  Pipes 
Newark,  DE  19711 

1  University  of  Denver 

Denver  Research  Institute 
ATTN:  Mr.  R.  F.  P.echt 
2390  S.  University  Blvd. 
Denver,  CO  80210 
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2  University  of  Florida 

Department  of  Engineering 
Sciences 

ATTN:  Dr.  R.  L.  Sierakowski 
Dr.  L.  E.  Malvern 
Gainesville,  FL  32601 

1  University  of  Oklahoma 
School  of  Aerospace, 
Mechanical  and  Nuclear 
Engineering 
ATTN:  Dr.  C.  W.  Bert 
Norman,  OK  7306V 

1  Commander 

Armament  R&D  Center 
US  Army  AMCCOM 
ATTN:  DRSMC-TDC(D) 

Dover,  NJ  07801 


Copies  Organization 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP ,  H.  Cohen 

Cdr ,  USATECOM 
ATTN:  DRSTE-TO-F 

Dir,  MTD 

ATTN:  Mr.  W.  Pleas 

Mr.  S.  Keithley 

Cdr,  CROC,  AMCCOM 
ATTN :  DRSMC-CLB-PA 
DRSMC-CLN 
PRSMC-CLJ-L 
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2  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Bennlng,  GA  31905 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-RBL 
Redstone  Arsenal,  AL  35898 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-YDL 
Redstone  Arsenal,  AL  35898 

3  Commander 

Naval  Sea  Systems  Command 
Navy  Department 
ATTN:  SEA-62R 
SEA-62Y 
SEA  9961 

Washington,  DC  20360 

1  HQ,  PACAF/DOOQ 
Hickam  AFB 
Honolulu,  HA  96853 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1.  BRL  Report  Number^ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


